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I. INTRODUCTION 

-xf ■ 

Since the observation of the current expansion of the Universe which has apparently accel- 
erated in the recent past, the anomalies found in the cosmic microwave background (CMB) and 
the large structures observations it becomes obvious that a pure Friedmann-Lemaitre-Robertson- 
q \ Walker (FLRW) cosmology should be amended. 

Bianchi type I cosmological models are the simplest anisotropic Universe models playing an 
important role in understanding essential features of the Universe. In this class of models it is 
possible to accommodate the presence of cosmic strings. In the last time the string cosmological 
models have been used in attempts to describe the early Universe and to investigate anisotropic 
dark energy component including a coupling between dark energy and a perfect fluid (dark matter) 
[Q]]. Cosmic strings are one dimensional topological defects associated with spontaneous symme- 
try breaking in gauge theories. Their presence in the early Universe can be justified in the frame 
of grand unified theories (GUT). 

The object of this paper is to investigate a Bianchi type I string cosmological model in the 
presence of a magnetic flux. The inclusion of the magnetic field is motivated by the observational 
cosmology and astrophysics indicating that many subsystems of the Universe possess magnetic 
fields (see e. g. the reviews [HLH] and references therein). 

In the following section we introduce a system of cosmic string and magnetic field in the 
Bianchi type I cosmology presenting some of its general features. In Section III we introduce a few 
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plausible assumptions usually accepted in the literature and some exact solutions are produced. In 
the last section we present some conclusions and perspectives. 



H. FUNDAMENTAL EQUATIONS AND GENERAL SOLUTIONS 

We consider the gravitational filed given by an anisotropic Bianchi type I (BI) metric 

ds 2 = al(dx°) 2 -a 2 1 (dx 1 ) 2 -4(dx 1 ) 2 -al(dx i ) 2 , (2.1) 

with ciq = 1 , x° = ct and c = 1 . The metric functions a,- (z = 1 , 2, 3) are the functions of time t only. 
Einstein's gravitational field equation for the BI space-time has the form 

02 + a 1 + 02a 1= ^ (22a) 

Cl2 Cl3 Cl2 

— + — + —— = kT 2 2 , (2.2b) 

$3 U\ CI3 U\ 

— + — + —— = kTJ, (2.2c) 

(l\ CI2 Cl\ 0-2 



a 1 a2 + a2a 1 + a 1 a 1= _ ^ 

CL\ 0-2 Cl2 a 3 a 3 a \ 

Here k is the Einstein gravitational constant and over-dot means differentiation with respect to 
t. The energy momentum tensor for a system of cosmic string and magnetic field in a comoving 
coordinate is given by 

where p is the rest energy density of strings with massive particles attached to them and can be 
expressed as p = p p + X, where p p is the rest energy density of the particles attached to the strings 
and X is the tension density of the system of strings HLBLIg], which may be positive or negative. 
Here Uj is the four velocity and Xj is the direction of the string, obeying the relation 

U(u l = —XiX* = 1 , iijx' = 0. (2.4) 

In (12.31) E^v is the electromagnetic field given by Lichnerowich |0] 



// \ 2 (u^u v --Sjj^j -h^tf 



(2.5) 



Here p is a constant characteristic of the medium and called the magnetic permeability. Typically 
p differs from unity only by a few parts in 10 5 (p > 1 for paramagnetic substances and p < 1 for 
diamagnetic). In (12.51) is the magnetic flux vector defined by 



1 

where *F„ V is the dual electromagnetic field tensor defined as 



h ll = —*F Vil u v ', (2.6) 
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Here F a $ is the electromagnetic field tensor and £^ V a/3 is tne totally anti- symmetric Levi-Civita 
tensor with £0123 = +1. Here the co-moving coordinates are taken to be u° = 1, u l = u 2 = u 3 = 
0. We choose the incident magnetic field to be in the direction of x-axis so that the magnetic 
flux vector has only one nontrivial component, namely h\ 7^ 0. In view of the aforementioned 
assumption from (12.61) one obtains F\ 2 = ^13 = 0. We also assume that the conductivity of the fluid 
is infinite. This leads to Fq\ = Fq 2 = F03 = 0. Thus We have only one non-vanishing component 
of Fjxv which is F23. Then from the first set of Maxwell equation 

F liv,p+ F vP;^+ F Pfi;v =° 5 (2-8) 
where the semicolon stands for covariant derivative, one finds 



F 23 = J, J = const. (2.9) 
Then from (12.61) in account of (12.71) one finds 

h l = ^^. (2.10) 

M a 2#3 

Finally, for one finds the following non-trivial components 



Taking the string along x l direction and using co-moving coordinates we have the following 
components of energy momentum tensor ||8|] : 

<j2 2 

T} = A + —4, (2.12b) 



J 2 



2jU T 2 ' 
2/1 T 2 



T 2 2 = (2.12c) 



a2 2 
^3 _ J "1 



where we used the definition 



n = <2.i2d) 



T = aia2^3- (2.13) 



It is indeed the volume scale of the BI space-time, i.e., x 
In view of T 2 2 = from (|2.2bl) . (I2.2cl) one finds 



a 2 = a 3 Dexp(xJ (2.14) 

with D and X being integration constants. Due to anisotropy of the source filed, in order to solve 
the remaining Einstein equation we have to impose some additional conditions. Here we give two 
different conditions. It can be shown that the metric functions can be expressed in terms of T. So 
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let us first derive the equation for x. Summation of Einstein Eqs. (I2.2al) . (I2.2bl) . (I2.2cl) and 3 times 
(l2T2dl) gives 

T 2 V 3jU T 



± = i K (p + A + ±l^). (2.15) 



Let us demand the energy-momentum to be conserved, i.e., T^. v = 0, which in our case takes the 
form 

--(tT °) - *T} - ^T 2 2 - ^r 3 3 = 0. (2.16) 
After a little manipulation from (|2.16l) one obtains 

p + -p-— A=0. (2.17) 

III. SOME EXAMPLES AND EXPLICIT SOLUTIONS 

In the literature there exists a number of relations between p and A , the simplest one being a 
proportionality relation: 

p = al (3.1) 
with the most usual choices of the constant a 

1 geometric string 

a = { 1 + co CO > 0, p string or Takabayasi string (3.2) 
-1 Reddy string. 

As it was mentioned earlier, to solve the Einstein equations completely, we need to impose 
some additional conditions. Here we consider the following two. 



A. Case 1 

This condition was first used by Bali IllOn . Following him let us assume that the expansion 
scalar (0) in the model is proportional to the eigenvalue o\ of the shear tensor o^. For the BI 
space-time we have 

8 = ^1 + ^ + ^1^ (3.3) 

Cl\ Q.2 CI3 

<& = - 5 (4- + - + -). (3.4) 
Writing the aforementioned condition as 

e=3nol (3.5) 

one comes to the following relation 

a\ = Z{a2a-i) , (3.6) 
where N = — (n+l)/ {An + 1) being the proportionality constant and Z is the integration constant. 
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From (|2.14l) and (13.61) after some manipulation for the metric functions one finds [H] 



a i 
"3 



z \/(N+l) T N/(N+l) 
T\1/2(J\H-1) 



In this case Eq. (12.171) takes the form 



D(-)' ' 'rap 

1 /T N l/2(AT+l) 

(-) ; exp 



■x 




df 




.2 


I 


x . 


j 




X 


f df 




2 


J ' 


T - 



(3.7a) 
(3.7b) 

(3.7c) 



N 



(3.8) 



Eq. (12.151) now reads 



t = ^K( P + X)T + X^ N - 1 ^ N+1 \ where X = K ^Z 2 /^ +1 ) 



(3.9) 



We will see later, the right hand side of (13.91 ) is the function of T, hence can be written as 

z = 3 r (x). (3.10) 
Eq. (13.101) admits first integral which can be written as 

t= y/2[E-U{x)\ 

with 



U(x) = -J ?(x)dT. 



(3.11) 
(3.12) 



The expression (13.121) can be viewed as potential, while £ as energy level. A detailed analysis of 
this mechanism can be seen in, e.g., Ill ill . 

Let us now study Eqs. (13.81 ) and (13.91 ) for different equations of state. Assuming the relation 
(13.11) between the pressure of the perfect fluid p and the tension density A from Eq. (13.81) one finds 



( 



N 



with the solution 

while the equation for % reads 



p =C T< N 



+1) 



3 1 N N-l 



(3.13) 



(3.14) 



(3.15) 



This equation can be set in the following form 



• 3KCo( a + l)(N+l) /a( u XJiV+l) 

N+a(N+\) N 



(3.16) 



where £q is the integration constant and is related to £ as £q = 2£. 
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FIG. 1 : Potential corresponding to the 
different equations of state in absence 
of a magnetic field. 




FIG. 2: Evolution of the Universe for 
different equations of state in absence 
of a magnetic field. 




FIG. 3: Potential corresponding to the 
different equations of state in pres- 
ence of a magnetic field. 



FIG. 4: Evolution of the Universe for 
different equations of state in pres- 
ence of a magnetic field. 



In Figs. [Hand [3] we have illustrated the potential corresponding to the different equations of 
state without and with the magnetic field, respectively. The Figs. [2] and |4] show the evolution of % 
for different cases. Here "G", "P" and "R" stand for geometric string, p string and Reddy string, 
respectively. The reason to illustrate figures with and without magnetic field is to show the role 
of magnetic field. As one sees, in all cases T might be zero at the initial stage of evolution, thus 
giving rise to the initial singularity in one hand, % is not bound from above which means in all 
three cases we have ever expanding Universe. But introduction of magnetic field into the system 
results in rapid growth of x. In numerical analysis we used the following value for the problem 
parameters: JC=l,ft)=l,./V = 4, Z=l.In case of string only we set J = 0, otherwise 3=1. For 
magnetic permeability we choose fx = 1.00001 and p, = 0.99999, respectively. Since it doesn't 
make any significant change in the behavior of T, we illustrate only case with p, = 1.00001. The 
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initial value for x is taken to be t(0) = 0.0001 and corresponding first derivative is calculated from 
(I3~i6l) at £n = 1. 



B. Case 2 

The second case was proposed by Barrow II 1211 . Here we first introduce the generalized Hubble 
parameter 

1 ,d\ a,2 (hx It ,„ „„ N 

H = - (— + - + - ) = --, (3.17) 
3 a\ ci2 ay 3 X 

and two relative shear anisotropy parameters by 

1 ,a,\ Qi \ 

R = -{---), (3.18a) 
H a\ a,2 

1 ,a,\ d?.\ 
ti a\ «3 

When R = S = the universe will be the isotropic flat Friedman universe. If one sets R = const., 
then in view of (13.171) one finds, 

fl!=D a2^ /3 . (3.19) 
Then, in view of (|2.13l) . (12.141) we find the following expressions for metric functions 

ai = {EftDfftTp+WftJSZ, (3.20a) 

a 2 = (zV£) 1/3 -r (3 -* )/9 e f/ S (3.20b) 

a 3 = (D D 2 )- l / 3 TC- R y 9 e=¥l*. (3.20c) 
In this case Eq. (12.171) takes the form 

P+(P-^A)1-|X=0, (3.21) 
whereas, for x in this case we have 

- = lK(p + X + ^(D 2 D) 2 ^x^- i5 y 9 e^fT). (3.22) 
x 2 V 3/1 / 

Comparing this equation with the corresponding one (|3.15l) from Case 1, it results that in Case 
2 the equation of evolution for x is more intricate. Indeed Eq. (|3.22l) involves the function x in an 
integral and needs further study. A detailed analysis is presently under investigation and will be 
reported in a forthcoming paper lfl3ll . 

In the present paper we confine ourselves to a preliminary task considering this equation in an 
asymptotic regime, for large t. To be aware of possible asymptotic behavior of p and x for large t, 
let us assume again relation (13.11) and Eq. (|3.21l) becomes 

I, ,3.23) 
p v 9a 'x 3a x 

which shows that both sides of this equation are the same function of time. As in the Case 1 Eq. 
(13.231) allows a behavior of the pressure of the perfect fluid p decreasing in time as an inverse 
power of t: 

p=C l t-P, (3.24) 
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C\ and /3 > being some constants. Correspondingly, for t we get the equation 
with the solution 

3X _9ap_ 

f + C 2 f~3+2K-9« . (3.26) 



9ap+3 + 2R-9ct 

It is interesting to observe that Eq. (13.231 ) allows even an exponential behavior of p for large V. 

p = C^e-V (3.27) 

which means that the both sides of Eq. (13.231) are a negative constant (—7). The corresponding 
behavior of % from Eq. (|3.23l) is 

X 9ay 

+ C 4 e~3+2*-9a' (3.28) 



3ay 



with C4 another arbitrary constant. 

The full study of Eq. (13.221) . which contains the parameters characterizing the cosmic string 
and magnetic field, for large t using the asymptotic behaviors (13.241) . (13.251) and (13.271) . (13.281) will 
be presented elsewhere II 1 311 . 



IV. CONCLUSIONS 



In the present paper we investigated in the frame of Bianchi type I models a string cosmological 
model in the presence of a magnetic field. We used some tractable assumptions concerning the 
parameters entering the model. 

In the case of a proportionality between the trace of the expansion tensor 9 and the eigenvalue 
o\ of the shear tensor we are able to get explicit analytic solutions. Setting the relative anisotropy 
parameter R = const the model is more involved, but sufficiently interesting to deserve further 
study. 



Acknowledgments 



This work is supported in part by a joint Romanian-LIT, JINR, Dubna Research Project, theme 
no. 05-6-1060-2005/2010. 



[1] Koivisto T. and Mota D.F., Phys. Lett. B, 644, 104 (2007); Koivisto T. and Mota D.F. , Phys. Rev. D, 

75, 023518 (2007). 

[2] Grasso D. and Rubinstein H.R., Phys. Rep., 348, 163 (2001). 

[3] Widrow L.M., Rev. Mod. Phys., 74, 775 (2002). 

[4] Letelier, PS., Phys. Rev. D, 28, 1424 (1983). 

[5] Pradhan A., Yadav A.K., Singh R.P and Singh V.K., Astrophys. Space ScL, 312, 145 (2007). 

[6] Khadekar G.S., Tade S.D., Astrophys. Space Sci. 310, 47 (2007)]. 



String cosmological model in the presence of a magnetic flux 



9 



[7] Lichnerowicz A., Relativistic Hydrodynamics and Magnetohydrodynamics, (Benjamin, New York, 
1967). 

[8] Saha Bijan, J. Astrophys. Space Sci. 2005. 299 No. 1. 1149-1581 (2005). 

[9] Saha Bijan, Phys. Rev. D. 64H2350LJ(2001). 

[10] Bali R. Int. J. Theor. Phys. 25, 755 (1986). 

[11] Bijan Saha and Todor Boyadjiev, Phys. Rev. D 69. 11240101 (2004) [arXiv: |gFqc70311045| . 

[12] Barrow J.D., Phys. Rev. D. 55 7451 (1997). 

[13] Saha B. and Visinescu M., in preparation. 



